Phonon-dressed Mollow triplet in the regime of cavity-QED 
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We study the resonance fluorescence spectra of a driven quantum dot placed inside a high Q 
semiconductor cavity and interacting with an acoustic phonon bath. The dynamics is calculated 
using a time-convolutionless master equation obtained in the polaron frame. We demonstrate pro- 
nounced spectral broadening of the Mollow sidebands through cavity-emission which, for small 
cavity-coupling rates, increases quadratically with the Rabi frequency. However, for larger cavity 
coupling rates, this broadening dependence is found to be more complex. This field-dependent 
Mollow triplet broadening is primarily a consequence of the triplet peaks sampling different parts 
of the asymmetric phonon bath, and agrees directly with recent experiments with semiconductor 
micropillars. The influence from the detuned cavity photon bath and multi-photon effects is shown 
to play a qualitatively important role on the fluorescence spectra. 

PACS numbers: 42.50.Ct, 78.67.Hc, 78.55.-m 



Recent developments in semiconductor cavity systems 
exploit several fields of research, including semiconductor 
optics, quantum information processing and cavity-QED 
(quantum electrodynamics). Target quantum sources 
such as indistinguishable photons are important in tech- 
nological applications [ljj, and require robust physics de- 
signs and an understanding of how to manipulate dynam- 
ical processes at a very fundamental level. Unique quan- 
tum processes can be obtained when the cavity and quan- 
tum dot (QD) are near the strong coupling regime 0>I3, 
which introduces quantum mechanical features such as 
photon antibunching Recent experimental studies 
involving coherent excitation have focussed on resonance 
fluorescence of a QD coupled to a cavity mode [g- 7] . Also 
of interest is the study of off-resonant QD cavity systems 
that can be used to monitor resonance fluorescence of 
the emitter [§[. Very large optical dipole moments rela- 
tive to atomic systems and strong interaction effects arc 
the salient features of these QDs. Yet, the very nature 
of the semiconductor QDs can also introduce excitation- 
dependent dephasing and decoherence effects. 

In 1969, Mollow demonstrated that the fluorescence 
spectrum of a two-level atom has additional sidebands 
in addition to the central coherent peak when the atom 
is driven by a strong coherent single- mode field 0, |T(| ■ 
This effect can be attributed to the emergence of "dressed 
states" — an infinite series of doublets or Floquet states, 
whose positions and spectral widths are well known for 
large driving fields. From both a practical and fundamen- 
tal level, it is interesting to study such well known atomic 
optics effects when a semiconductor QD is strongly driven 
by a coherent laser field. For pulsed-excited systems, 
excitation-induced dephasing (EID) emerges due to the 
interaction of the laser-induced dressed states with the 
underlying acoustic phonon reservoir llHl3j. For con- 
tinuous (cw) excitation, it is not clear what the role 
of phonon coupling will be, as it is often argued that 
the dephasing processes are due to the "non-Markovian" 
processes, especially important for short pulses and fast 



timescales. Thus one may expect that cw excitation 
does not suffer much from phonon-induced damping, and 
indeed recent experimental works using QD resonance 
fluorescence [51421 suggest negligible influence from the 
phonon bath. Nevertheless, if a QD system is driven 
strongly enough, a phonon-induced dephasing process 



may occur 14[. In recent experiments [15j, cw reso- 



nance fluorescence emission from single QD placed in- 
side a high-quality micropillar cavity was studied. In 
contrast to all previous works on cw-excited semicon- 
ductor systems, it was found that increasing excitation 
power results in a Mollow triplet spectra with a system- 
atic spectral sideband broadening that is proportional to 
the square of the Rabi fre quen cy — which was interpreted 
as a signature of EID ll|, ll2[ . Since the QD was placed 
in the regime of cavity-QED, which enables significant 
cavity coupling, the origin of this Mollow sideband de- 
phasing is not at all clear. 

In this work, we develop a nonperturbative theory 
of cavity-QED coupling in a longitudinal acoustic (LA) 
phonon bath. We apply this theory to connect directly to 
the EID observations on the Mollow triplet, and we intro- 
duce the phenomenon of: phonon-dressed Mollow triplet 
in the regime of cavity-QED. We develop a model that de- 
scribes both electron-acoustic phonon coupling and cav- 
ity photon coupling to all orders. Specifically, we describe 
the process of electron-photon and electron-phonon cou- 
pling by using a polaron transform [lij , which, in the ap- 
propriate limits, formally recovers the independent boson 
model or the Jaynes-Cumming model. The polaron trans- 
form allows one to eliminate the QD phonon interaction 
part from the system Hamiltonian at the expense of mod- 
ified Rabi frequency, and the dot-cavity coupling that 
is now dressed by coherent phonon displacement opera- 
tors [16, 17 1 . We then perform time local projection oper- 
ations and derive a straightforward time-convolutionless 
master equation that treats the residual system-bath in- 
teraction to second order [l8| , while treating cavity pho- 
ton effects to all orders. To the best of our knowledge, 
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this is the first study of its kind to explore electron- LA- 
phonon interactions and cavity photon interactions to all 
orders. 

Working in a frame rotating with respect to the laser 
pump frequency, the model Hamiltonian is 

H = H(oj x — ujl)(t + (j~ + h(uj c — ui]f)a) a 
+ hg(a + a + aV~) + hfl(a + + <r~) 

+ o + a- hX k(b k + b{) + h"kb{b k , (1) 

k k 

where bk{b* k ) are the annihilation and creation operators 
of the phonon reservoir, a is the leaky cavity mode anni- 



hilation operator, a + ,a (and o~ z 



a o~ + ) arc 



the Pauli operators of the electron- hole pair ("exciton"), 
S a L = u a — lul (a = x,c) are the detunings of the ex- 
citon and cavity from the coherent pump laser, and g is 
the cavity-exciton coupling strength. We next perform 
a polaron-transformation that introduces a renormalized 
Rabi frequency and dot-cavity coupling strength 17 
H' = cxp(S*)iJexp(-5), with S = a+a~ £ fc ^{b\-b k 
The transformed Hamiltonian has three contributions, 



n sys 



h(5 xL - Ap)<7 + <7 + M cL a)a + (B)X g , (2a) 
H Lth = ^2^kb{h, 
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where X„ 



ih[g(a + a — a^ a ) 

b% c g = U B + 

Note that (B) = 



%(aV 



(2b) 
(2c) 
X u = 



-il(a- -*+)], B ± = exp[±£ fe ^(6 fe - 

B_ - 2(B)), and ( u = ±(B + - 
(B + ) = (B_) and the polaron shift is 



given by Ap = J2 k ^T- Without loss of generality, we as- 
sume that the polaron shift is implicitly included in tu x . 
It is important to note the slightly different definition of 
the system Hamiltonian in Eq. (f2a|) . The usual decom- 
position of the system Hamiltonian to include only the 
noninteracting QD and cavity parts is known to fail the 
detailed balance condition in general, which is caused by 
internal coupling effects [20| . The transformed system 
Hamiltonian above leads to the correct form of the den- 
sity operator and preserves detailed balance. 

Next, we introduce dissipative terms that describe the 
radiative decay of the QD exciton and the cavity mode. 
They are included as Liouvillian superoperators acting 
on the reduced system density matrix 21|: 



L{p) = -^-(2a po + — o + a p — pa + a ) 

r' 

+ n{2apa) — o) ap — pa) a) + — (a z po~ z — p), (3) 

where 2k is the cavity decay rate, T' is the pure dephas- 
ing rate, and T x = T X {B) 2 is the radiative decay rate. 
The latter decay has an additional rcnormalization by a 



factor of (B) 2 resulting from the polaron transform. The 
physical origin of this can be traced to the dephasing 
of the optical dipole moment resulting from the Franck- 
Condon displacement of the excited state 22J. In writ- 
ing the phonon induced renormalization in this form, we 
have assumed a separation of the timescales associated 
with the radiative and phonon processes which amounts 
to replacing the coherent phonon displacement operators 
B± by (B) = (B±), where the angular brackets denote 
trace over the thermal phonon reservoir. This can be jus- 
tified by noting that phonon dynamics is extremely fast 
and happens at picosecond timescales, i.e., much faster 
than the typical radiative processes which occur. 

We now derive a time-local master equation [l8| using 
the time-convolutionless approach for the cavity-QED 
system density operator in the interaction picture, de- 
scribed by H' sys , and to second order in exciton-photon- 
phonon coupling H' int . The phonon reservoir is consid- 
ered to be stationary and in thermal equilibrium and 
factorized initial conditions are assumed. The time- 
convolutionless master equation for the cavity-QED sys- 
tem density operator is then given by 



dp_]_ 

dt ih 



[H' sys ,p{t)]+L{p)-^ f dr ( 4 ) 

7n—g,u 

x {G m (T)[X m ,e- iH *v° T / n X m e iH »v* T / h p(t)]+H.c.), 



where G g / U (t) = (( g / u (t)( g /u{0))- These polaron Green 
functions are well known and are defined as 0, [ItJ 
G g (t) = (B) 2 (coshk6(i)] - 1), G u (t) = (B) 2 sinh[0(i)], 
with (f>(t) = / °° d(jJ : ^jr-{coth({3huj/2) cos(wi) — isin(art)]. 
The (B) and the polaron-shift Ap can now be written 
as (B) = cxp(-i/ °°dcj^coth(^?iu;/2)) and A P = 

f n °° duj^^-. This model describes the electron-LA- 
phonon interaction via a deformation potential cou- 
pling which is the primary source of dephasing in self- 
assembled QDs [23 1 . The spectral function that con- 
siders electron-phonon interactions via a deformation 
potential coupling can now be written as J(u>) = 
a p uj 3 cxp(— uj 2 /2uj 2 ) where we use w& = 1 meV and 
ap/(27r) 2 = 0.06 ps 2 as typical numbers for InAs/GaAs 



QDs [24j, |25|. With these parameters, then the phonon- 
induced renormalization contribution, (B) = 0.84, at 
T = 10 K. With the polaron master equation, we 
can readily calculate the cavity-mode incoherent emis- 
sion, from: S ca v{u) oc limt^.ooRe{/ dt[(a(t + r)a^(t)) — 
(a(t + t)) (a t (i))]e i (" £ - -")*}. The master equation is nu- 
merically solved in a basis that can be truncated at any 
arbitrary photon and exciton state, which enables us to 
compute the weak excitation approximation results (one 
quanta limit) and the regime of multi-quanta cavity- 
QED. Since we consider cw excitation, and the phonon 
dephasing timescales are very fast, we can take t — > 00 
in Eq. @; we have checked this approximation to be rig- 
orously valid which confirms that that phonon processes 
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FIG. 1: (Color online) Cavity-emitted fluorescence spectra 
for the one-phonon model (blue-dashed) and the full polaron 
theory (red-solid). The cw Rabi drive strength Q — 0.05 meV, 
and is resonant with the exciton (5 x l = 0). The parameters 
were as follows: g = 0.02 meV, V' = 2 n&V, T x = 1 HbV, 
2k = 0.1 meV and to c — uj x = 0.28 meV. 
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FIG. 2: (Color online) Cavity-emitted fluorescence spectra 
with multi-photon processes (red-solid) compared to one- 
photon processes (blue-dashed) in the quantum dot cavity 
system. We consider phonons to all orders, and two differ- 
ent dot-cavity coupling parameters (a) g = 0.02 meV and (b) 
g = 0.05 meV. The cw Rabi drive and material parameters as 
the same as those used in Fig. [1] 



occur at timcscalcs substantially faster than the system 
dynamics. The two-time correlation function is obtained 
from the quantum regression formula (Io| . 

We first study the need for both multi-phonons and 
multi-photons. To be consistent with the recent Mollow 
triplet experiments on semiconductor micropillars fl5| . 
we consider a detuned cavity mode with u> x — cu c w 
0.3 meV, a cw resonantly-excited exciton, and emitted 
spectra obtained from cavity-emission, S cav . The radia- 
tive decay rate T x = 1 fi eV and the residual pure de- 
phasing rate, V = 2_u eV; these are similar to those 
used by Ulrich et al. [15| , and between those measured in 
Ref. [l[ and Ref. @ . In Fig. [TJ we plot the cavity-emitted 
spectra for the one-phonon model and compare it with 
the full polaron theory (phonons to all orders) . The one- 
phonon correlation function is obtained by expanding the 
phonon correlation functions G giU (t), and (B), to lowest 
order in the dot-phonon coupling. At low temperatures, 
the one-phonon expansion of the polaronic phonon cor- 
relation function is sufficient to describe the dynamics. 



However, with increasing temperatures, the one-phonon 
model starts to overestimate the effect of phonon cou- 
pling with respect to the polaronic description. This 
trend is in consistent with the pulse-excitation analysis 
of McCutcheon and Nazir 181, for T > 30 K; however, 



what is particularly surprising in our cw-cavity case is 
that multi-phonon effects are important at T = 10 K. 

In Fig. [5]we study the role of multi-photon processes in 
the QD cavity dynamics in the presence of phonons, for 
different values of dot-cavity coupling parameters. We 
find that a one-photon-correlation approximation fails 
in the present excitation regime, though a two-photon- 
correlation basis was found to be sufficient for all the cal- 
culations that follow. This confirms that the theory must 
include multi-phonon and multi-photon effects, even for 
a detuned cavity system at T = 10 K. With photons and 
phonons included, we recognize two remarkable features 
in the cavity-emitted spectra. First, there is a strong 
asymmetry for the on-resonance Mollow triplet, which is 
primarily due to the off-resonant cavity. Second, there is 
a pronounced cavity feeding effect, where the cavity mode 
is excited via phonon interactions, which is a similar ef- 



fect to incoherently excited QD-cavity systems [21|, [24- 
[29j , but much stronger in this coherent excitation regime. 
The asymmetries of the Mollow triplet are attributed to 
the subtle interplay of the cavity-photon bath, phonon 
bath and pure dephasing. 

In Fig. [3l wc show the cavity-emitted spectrum of 
the QD exciton under resonant excitation conditions 
($xL — 0) as a function of cw excitation Rabi frequency, 
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FIG. 3: (Color online) Cavity-emitted fluorescence spectra of 
a quantum dot placed under cw excitation of the exciton (with 
fix =0). We vary the Rabi frequency from 0.01 — 0.05 meV, 
and we consider two different quantum dot-cavity coupling 
parameters (a) g — 15 /ieV and (b) g = 50 fieV. The blue- 
dashed plots correspond to the cases with no phonon cou- 
pling, while the red-solid line corresponds to the case where 
the quantum dot is coupled to a phonon reservoir at a tem- 
perature T = 10 K. The other parameters are as in Fig. [TJ 
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excellent agreement with the very recent experiments of 
Ulrich et al. 
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FIG. 4: (Color online) FWHM of the Mollow sidebands as a 
function of square of the Rabi frequency (fi 2 ) extracted by a 
numerical fit from Fig. [3] in the presence of phonon coupling 
(red-solid) and with no phonon coupling (blue-dashed). The 
linear fit is a signature of excitation-induced dephasing due 
to the underlying phonon coupling to the quantum dot. 



for g = 15 /xeV and g = 50 /xeV. The spectra are nor- 
malized to unity, and, for clarity, the spectra have been 
shifted vertically as a function of Rabi field strength. We 
show spectra with and without phonon coupling, where 
the former assumes a phonon reservoir in thermal equi- 
librium at T = 10 K. As expected, we clearly find in- 
creasing splitting of the Mollow triplet with increasing 
Rabi frequency f2. We also note that the location of the 
sidebands changes due to the renormalization of the Rabi 
frequency in the presence of phonon coupling which re- 
duces the effective Rabi frequency; the Mollow splitting 
is slightly reduced as a result and the triplet resonances 
are significantly broader. In Fig. |4l we plot the corre- 
sponding numerically-extracted FWHM (full-width half- 
maximum) values as a function of the square of the Rabi 
frequency. We find a very good linear fit to the data 
which suggests that Mollow triplets are broadened with 
increasing f2, in agreement with recent experiments fl5T ]. 
However, for the case of g = 50 neV, the dependence is 
likely nonlinear for low values of f2, and in general de- 
phasing results from the interplay of both phonon bath 
coupling and cavity photon bath coupling. In addition, 
we obtain rich asymmetries of the Mollow triplet and 
substantial cavity feeding - two pronounced effects that 
are unique to the QD-cavity system. 

In summary, we have presented a theoretical and nu- 
merical study of the cavity-emitted resonance fluores- 
cence of a QD placed inside a high Q cavity. We de- 
veloped a polaron master equation formalism based on 
the time-convolutionlcss approach which enabled us to 
study the dynamical properties of the reduced density 
matrix of the QD-cavity system in response to both an 
external cw laser field and a bath of thermalized LA 
phonons. We applied this theory to study the influence 
of the Rabi frequency on cavity-emitted spectra and find 
strong numerical evidence of excitation-induced dephas- 
ing. Both multi-phonon and multi-photon effects arc 
shown to be important and our general findings are in 
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